Towards identification of postharvest fruit quality transcriptomic markers in Malus domestica.

Gene expression is highly impacted by the environment and can be reflective of past events that affected developmental processes. It is therefore expected that gene expression can serve as a signal of a current or future phenotypic traits. In this paper we identify sets of genes, which we call Prognostic Transcriptomic Biomarkers (PTBs), that can predict firmness in Malus domestica (apple) fruits. In apples, all individuals of a cultivar are clones, and differences in fruit quality are due to the environment. The apples transcriptome responds to these differences in environment, which makes PTBs an attractive predictor of future fruit quality. PTBs have the potential to enhance supply chain efficiency, reduce crop loss, and provide higher and more consistent quality for consumers. However, several questions must be addressed. In this paper we answer the question of which of two common modeling approaches, Random Forest or ElasticNet, outperforms the other. We answer if PTBs with few genes are efficient at predicting traits. This is important because we need few genes to perform qPCR, and we answer the question if qPCR is a cost-effective assay as input for PTBs modeled using high-throughput RNA-seq. To do this, we conducted a pilot study using fruit texture in the 'Gala' variety of apples across several postharvest storage regiments. Fruit texture in 'Gala' apples is highly controllable by post-harvest treatments and is therefore a good candidate to explore the use of PTBs. We find that the RandomForest model is more consistent than an ElasticNet model and is predictive of firmness (r2 = 0.78) with as few as 15 genes. We also show that qPCR is reasonably consistent with RNA-seq in a follow up experiment. Results are promising for PTBs, yet more work is needed to ensure that PTBs are robust across various environmental conditions and storage treatments.

A chromosome-scale assembly for 'd'Anjou' pear.

Cultivated pear consists of several Pyrus species with Pyrus communis (European pear) representing a large fraction of worldwide production. As a relatively recently domesticated crop and perennial tree, pear can benefit from genome-assisted breeding. Additionally, comparative genomics within Rosaceae promises greater understanding of evolution within this economically important family. Here, we generate a fully phased chromosome-scale genome assembly of P. communis 'd'Anjou.' Using PacBio HiFi and Dovetail Omni-C reads, the genome is resolved into the expected 17 chromosomes, with each haplotype totaling nearly 540 Megabases and a contig N50 of nearly 14 Mb. Both haplotypes are highly syntenic to each other and to the Malus domestica 'Honeycrisp' apple genome. Nearly 45,000 genes were annotated in each haplotype, over 90% of which have direct RNA-seq expression evidence. We detect signatures of the known whole-genome duplication shared between apple and pear, and we estimate 57% of d'Anjou genes are retained in duplicate derived from this event. This genome highlights the value of generating phased diploid assemblies for recovering the full allelic complement in highly heterozygous crop species.

Transcriptomic approach to uncover dynamic events in the development of mid-season sunburn in apple fruit.

Apples grown in high heat, high light, and low humidity environments are at risk for sun injury disorders like sunburn and associated crop losses. Understanding the physiological and molecular mechanisms underlying sunburn will support improvement of mitigation strategies and breeding for more resilient varieties. Numerous studies have highlighted key biochemical processes involved in sun injury, such as the phenylpropanoid and reactive oxygen species (ROS) pathways, demonstrating both enzyme activities and expression of related genes in response to sunburn conditions. Most previous studies have focused on at-harvest activity of a small number of genes in response to heat stress. Thus, it remains unclear how stress events earlier in the season affect physiology and gene expression. Here, we applied heat stress to mid-season apples in the field and collected tissue along a time course-24, 48, and 72 h following a heat stimulus-to investigate dynamic gene expression changes using a transcriptomic lens. We found a relatively small number of differentially expressed genes (DEGs) and enriched functional terms in response to heat treatments. Only a few of these belonged to pathways previously described to be involved in sunburn, such as the AsA-GSH pathway, while most DEGs had not yet been implicated in sunburn or heat stress in pome fruit.

Building a foundation for gene family analysis in Rosaceae genomes with a novel workflow: A case study in Pyrus architecture genes.

The rapid development of sequencing technologies has led to a deeper understanding of plant genomes. However, direct experimental evidence connecting genes to important agronomic traits is still lacking in most non-model plants. For instance, the genetic mechanisms underlying plant architecture are poorly understood in pome fruit trees, creating a major hurdle in developing new cultivars with desirable architecture, such as dwarfing rootstocks in European pear (Pyrus communis). An efficient way to identify genetic factors for important traits in non-model organisms can be to transfer knowledge across genomes. However, major obstacles exist, including complex evolutionary histories and variable quality and content of publicly available plant genomes. As researchers aim to link genes to traits of interest, these challenges can impede the transfer of experimental evidence across plant species, namely in the curation of high-quality, high-confidence gene models in an evolutionary context. Here we present a workflow using a collection of bioinformatic tools for the curation of deeply conserved gene families of interest across plant genomes. To study gene families involved in tree architecture in European pear and other rosaceous species, we used our workflow, plus a draft genome assembly and high-quality annotation of a second P. communis cultivar, 'd'Anjou.' Our comparative gene family approach revealed significant issues with the most recent 'Bartlett' genome - primarily thousands of missing genes due to methodological bias. After correcting assembly errors on a global scale in the 'Bartlett' genome, we used our workflow for targeted improvement of our genes of interest in both P. communis genomes, thus laying the groundwork for future functional studies in pear tree architecture. Further, our global gene family classification of 15 genomes across 6 genera provides a valuable and previously unavailable resource for the Rosaceae research community. With it, orthologs and other gene family members can be easily identified across any of the classified genomes. Importantly, our workflow can be easily adopted for any other plant genomes and gene families of interest.

GEMmaker: process massive RNA-seq datasets on heterogeneous computational infrastructure.

BACKGROUND: Quantification of gene expression from RNA-seq data is a prerequisite for transcriptome analysis such as differential gene expression analysis and gene co-expression network construction. Individual RNA-seq experiments are larger and combining multiple experiments from sequence repositories can result in datasets with thousands of samples. Processing hundreds to thousands of RNA-seq data can result in challenges related to data management, access to sufficient computational resources, navigation of high-performance computing (HPC) systems, installation of required software dependencies, and reproducibility. Processing of larger and deeper RNA-seq experiments will become more common as sequencing technology matures. RESULTS: GEMmaker, is a nf-core compliant, Nextflow workflow, that quantifies gene expression from small to massive RNA-seq datasets. GEMmaker ensures results are highly reproducible through the use of versioned containerized software that can be executed on a single workstation, institutional compute cluster, Kubernetes platform or the cloud. GEMmaker supports popular alignment and quantification tools providing results in raw and normalized formats. GEMmaker is unique in that it can scale to process thousands of local or remote stored samples without exceeding available data storage. CONCLUSIONS: Workflows that quantify gene expression are not new, and many already address issues of portability, reusability, and scale in terms of access to CPUs. GEMmaker provides these benefits and adds the ability to scale despite low data storage infrastructure. This allows users to process hundreds to thousands of RNA-seq samples even when data storage resources are limited. GEMmaker is freely available and fully documented with step-by-step setup and execution instructions.

PlantTribes2: Tools for comparative gene family analysis in plant genomics.

Plant genome-scale resources are being generated at an increasing rate as sequencing technologies continue to improve and raw data costs continue to fall; however, the cost of downstream analyses remains large. This has resulted in a considerable range of genome assembly and annotation qualities across plant genomes due to their varying sizes, complexity, and the technology used for the assembly and annotation. To effectively work across genomes, researchers increasingly rely on comparative genomic approaches that integrate across plant community resources and data types. Such efforts have aided the genome annotation process and yielded novel insights into the evolutionary history of genomes and gene families, including complex non-model organisms. The essential tools to achieve these insights rely on gene family analysis at a genome-scale, but they are not well integrated for rapid analysis of new data, and the learning curve can be steep. Here we present PlantTribes2, a scalable, easily accessible, highly customizable, and broadly applicable gene family analysis framework with multiple entry points including user provided data. It uses objective classifications of annotated protein sequences from existing, high-quality plant genomes for comparative and evolutionary studies. PlantTribes2 can improve transcript models and then sort them, either genome-scale annotations or individual gene coding sequences, into pre-computed orthologous gene family clusters with rich functional annotation information. Then, for gene families of interest, PlantTribes2 performs downstream analyses and customizable visualizations including, (1) multiple sequence alignment, (2) gene family phylogeny, (3) estimation of synonymous and non-synonymous substitution rates among homologous sequences, and (4) inference of large-scale duplication events. We give examples of PlantTribes2 applications in functional genomic studies of economically important plant families, namely transcriptomics in the weedy Orobanchaceae and a core orthogroup analysis (CROG) in Rosaceae. PlantTribes2 is freely available for use within the main public Galaxy instance and can be downloaded from GitHub or Bioconda. Importantly, PlantTribes2 can be readily adapted for use with genomic and transcriptomic data from any kind of organism.

A phased, chromosome-scale genome of 'Honeycrisp' apple (Malus domestica).

The apple cultivar 'Honeycrisp' has superior fruit quality traits, cold hardiness, and disease resistance, making it a popular breeding parent. However, it suffers from several physiological disorders, production, and postharvest issues. Despite several available apple genome sequences, understanding of the genetic mechanisms underlying cultivar-specific traits remains lacking. Here, we present a highly contiguous, fully phased, chromosome-level genome of 'Honeycrisp' apples, using PacBio HiFi, Omni-C, and Illumina sequencing platforms, with two assembled haplomes of 674 Mbp and 660 Mbp, and contig N50 values of 32.8 Mbp and 31.6 Mbp, respectively. Overall, 47,563 and 48,655 protein-coding genes were annotated from each haplome, capturing 96.8-97.4% complete BUSCOs in the eudicot database. Gene family analysis reveals most 'Honeycrisp' genes are assigned into orthogroups shared with other genomes, with 121 'Honeycrisp'-specific orthogroups. This resource is valuable for understanding the genetic basis of important traits in apples and related Rosaceae species to enhance breeding efforts.

Transcriptomics of Differential Ripening in 'd'Anjou' Pear (Pyrus communis L.).

Estimating maturity in pome fruits is a critical task that directs virtually all postharvest supply chain decisions. This is especially important for European pear (Pyrus communis) cultivars because losses due to spoilage and senescence must be minimized while ensuring proper ripening capacity is achieved (in part by satisfying a fruit chilling requirement). Reliable methods are lacking for accurate estimation of pear fruit maturity, and because ripening is maturity dependent it makes predicting ripening capacity a challenge. In this study of the European pear cultivar 'd'Anjou', we sorted fruit at harvest based upon on-tree fruit position to build contrasts of maturity. Our sorting scheme showed clear contrasts of maturity between canopy positions, yet there was substantial overlap in the distribution of values for the index of absorbance difference (I AD ), a non-destructive spectroscopic measurement that has been used as a proxy for pome fruit maturity. This presented an opportunity to explore a contrast of maturity that was more subtle than I AD could differentiate, and thus guided our subsequent transcriptome analysis of tissue samples taken at harvest and during storage. Using a novel approach that tests for condition-specific differences of co-expressed genes, we discovered genes with a phased character that mirrored our sorting scheme. The expression patterns of these genes are associated with fruit quality and ripening differences across the experiment. Functional profiles of these co-expressed genes are concordant with previous findings, and also offer new clues, and thus hypotheses, about genes involved in pear fruit quality, maturity, and ripening. This work may lead to new tools for enhanced postharvest management based on activity of gene co-expression modules, rather than individual genes. Further, our results indicate that modules may have utility within specific windows of time during postharvest management of 'd'Anjou' pear.

Comparative Analysis of the Apple Root Transcriptome as Affected by Rootstock Genotype and Brassicaceae Seed Meal Soil Amendment: Implications for Plant Health.

Brassicaceae seed meal (SM) soil amendment has been utilized as an effective strategy to control the biological complex of organisms, which includes oomycetes, fungi, and parasitic nematodes, that incites the phenomenon termed apple replant disease. Soil-borne disease control attained in response to Brassicaceae SM amendment is reliant on multiple chemical and biological attributes, including specific SM-generated modifications to the soil/rhizosphere microbiome. In this study, we conducted a comparative analyses of apple root gene expression as influenced by rootstock genotype combined with a seed meal (SM) soil amendment. Apple replant disease (ARD) susceptible (M.26) and tolerant (G.210) rootstocks cultivated in SM-amended soil exhibited differential gene expression relative to corresponding non-treated control (NTC) orchard soil. The temporal dynamics of gene expression indicated that the SM-amended soil system altered the trajectory of the root transcriptome in a genotype-specific manner. In both genotypes, the expression of genes related to plant defense and hormone signaling were altered in SM-amended soil, suggesting SM-responsive phytohormone regulation. Altered gene expression was temporally associated with changes in rhizosphere microbiome density and composition in the SM-treated soil. Gene expression analysis across the two rootstocks cultivated in the pathogen-infested NTC soil showed genotype-specific responses indicative of different defensive strategies. These results are consistent with previously described resistance mechanisms of ARD "tolerant" rootstock cultivars and also add to our understanding of the multiple mechanisms by which SM soil amendment and the resulting rhizosphere microbiome affect apple rootstock physiology. Future studies which assess transcriptomic and metagenomic data in parallel will be important for illuminating important connections between specific rhizosphere microbiota, gene-regulation, and plant health.

SHR4z, a novel decoy effector from the haustorium of the parasitic weed Striga gesnerioides, suppresses host plant immunity.

Cowpea (Vigna unguiculata) cultivar B301 is resistant to races SG4 and SG3 of the root parasitic weed Striga gesnerioides, developing a hypersensitive response (HR) at the site of parasite attachment. By contrast, race SG4z overcomes B301 resistance and successfully parasitises the plant. Comparative transcriptomics and in silico analysis identified a small secreted effector protein dubbed Suppressor of Host Resistance 4z (SHR4z) in the SG4z haustorium that upon transfer to the host roots causes a loss of host immunity (i.e. decreased HR and increased parasite growth). SHR4z has significant homology to the short leucine-rich repeat (LRR) domain of SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK) family proteins and functions by binding to VuPOB1, a host BTB-BACK domain-containing ubiquitin E3 ligase homologue, leading to its rapid turnover. VuPOB1 is shown to be a positive regulator of HR since silencing of VuPOB1 expression in transgenic B301 roots lowers the frequency of HR and increases the levels of successful SG4 parasitism and overexpression decreases parasitism by SG4z. These findings provide new insights into how parasitic weeds overcome host defences and could potentially contribute to the development of novel strategies for controlling Striga and other parasitic weeds thereby enhancing crop productivity and food security globally.

Temporal Dynamics of the Soil Metabolome and Microbiome During Simulated Anaerobic Soil Disinfestation.

Significant interest exists in engineering the soil microbiome to attain suppression of soil-borne plant diseases. Anaerobic soil disinfestation (ASD) has potential as a biologically regulated disease control method; however, the role of specific metabolites and microbial community dynamics contributing to ASD mediated disease control is mostly uncharacterized. Understanding the trajectory of co-evolutionary processes leading to syntrophic generation of functional metabolites during ASD is a necessary prelude to the predictive utilization of this disease management approach. Consequently, metabolic and microbial community profiling were used to generate highly dimensional datasets and network analysis to identify sequential transformations through aerobic, facultatively anaerobic, and anaerobic soil phases of the ASD process and distinct groups of metabolites and microorganisms linked with those stages. Transient alterations in abundance of specific microbial groups, not consistently accounted for in previous studies of the ASD process, were documented in this time-course study. Such events initially were associated with increases and subsequent diminution in highly labile metabolites conferred by the carbon input. Proliferation and dynamic compositional changes in the Firmicutes community continued throughout the anaerobic phase and was linked to temporal changes in metabolite abundance including accumulation of small chain organic acids, methyl sulfide compounds, hydrocarbons, and p-cresol with antimicrobial properties. Novel potential modes of disease control during ASD were identified and the importance of the amendment and "community metabolism" for temporally supplying specific classes of labile compounds were revealed.

Genome Sequence of Striga asiatica Provides Insight into the Evolution of Plant Parasitism.

Parasitic plants in the genus Striga, commonly known as witchweeds, cause major crop losses in sub-Saharan Africa and pose a threat to agriculture worldwide. An understanding of Striga parasite biology, which could lead to agricultural solutions, has been hampered by the lack of genome information. Here, we report the draft genome sequence of Striga asiatica with 34,577 predicted protein-coding genes, which reflects gene family contractions and expansions that are consistent with a three-phase model of parasitic plant genome evolution. Striga seeds germinate in response to host-derived strigolactones (SLs) and then develop a specialized penetration structure, the haustorium, to invade the host root. A family of SL receptors has undergone a striking expansion, suggesting a molecular basis for the evolution of broad host range among Striga spp. We found that genes involved in lateral root development in non-parasitic model species are coordinately induced during haustorium development in Striga, suggesting a pathway that was partly co-opted during the evolution of the haustorium. In addition, we found evidence for horizontal transfer of host genes as well as retrotransposons, indicating gene flow to S. asiatica from hosts. Our results provide valuable insights into the evolution of parasitism and a key resource for the future development of Striga control strategies.

Risk versus reward: host dependent parasite mortality rates and phenotypes in the facultative generalist Triphysaria versicolor.

BACKGROUND: Parasitic plants engage in a complex molecular dialog with potential host plants to identify a host and overcome host defenses to initiate development of the parasitic feeding organ, the haustorium, invade host tissues, and withdraw water and nutrients. While one of two critical signaling events in the parasitic plant life cycle (germination via stimulant chemicals) has been relatively well-studied, the signaling event that triggers haustorium formation remains elusive. Elucidation of this poorly understood molecular dialogue will shed light on plant-plant communication, parasitic plant physiology, and the evolution of parasitism in plants. RESULTS: Here we present an experimental framework that develops easily quantifiable contrasts for the facultative generalist parasitic plant, Triphysaria, as it feeds across a broad range of diverse flowering plants. The contrasts, including variable parasite growth form and mortality when grown with different hosts, suggest a dynamic and host-dependent molecular dialogue between the parasite and host. Finally, by comparing transcriptome datasets from attached versus unattached parasites we gain insight into some of the physiological processes that are altered during parasitic behavior including shifts in photosynthesis-related and stress response genes. CONCLUSIONS: This work sheds light on Triphysaria's parasitic life habit and is an important step towards understanding the mechanisms of haustorium initiation factor perception, a unique form of plant-plant communication.

"Stealth dissemination" of macrophage-tumor cell fusions cultured from blood of patients with pancreatic ductal adenocarcinoma.

Here we describe isolation and characterization of macrophage-tumor cell fusions (MTFs) from the blood of pancreatic ductal adenocarcinoma (PDAC) patients. The MTFs were generally aneuploidy, and immunophenotypic characterizations showed that the MTFs express markers characteristic of PDAC and stem cells, as well as M2-polarized macrophages. Single cell RNASeq analyses showed that the MTFs express many transcripts implicated in cancer progression, LINE1 retrotransposons, and very high levels of several long non-coding transcripts involved in metastasis (such as MALAT1). When cultured MTFs were transplanted orthotopically into mouse pancreas, they grew as obvious well-differentiated islands of cells, but they also disseminated widely throughout multiple tissues in "stealth" fashion. They were found distributed throughout multiple organs at 4, 8, or 12 weeks after transplantation (including liver, spleen, lung), occurring as single cells or small groups of cells, without formation of obvious tumors or any apparent progression over the 4 to 12 week period. We suggest that MTFs form continually during PDAC development, and that they disseminate early in cancer progression, forming "niches" at distant sites for subsequent colonization by metastasis-initiating cells.

A practical examination of RNA isolation methods for European pear (Pyrus communis).

OBJECTIVE: With the goal of identifying fast, reliable, and broadly applicable RNA isolation methods in European pear fruit for downstream transcriptome analysis, we evaluated several commercially available kit-based RNA isolation methods, plus our modified version of a published cetyl trimethyl ammonium bromide (CTAB)-based method. RESULTS: We confirmed previous work indicating chaotropic agent-based kits produced sufficient, high-quality RNA in freshly harvested, mature 'Bartlett' fruit. However, RNA isolation from 'd'Anjou' pear peel and especially cortical tissues of fruit stored for 11 months proved challenging to all but our modified CTAB-based method. Generally, more RNA was recovered from peel tissues than cortical tissues. Less toxic dithiothreitol was confirmed to be an acceptable reducing agent as its substitution for 2-mercaptoethanol often yielded high quality RNA. Finally, we present evidence that erroneous signals in the 5S region of Bioanalyzer RNA size plot histograms, that interfered with RNA integrity number calculation, were small RNA fragments that are reduced by simple cleanup procedures, not artifacts as previously reported.

Horizontal gene transfer is more frequent with increased heterotrophy and contributes to parasite adaptation.

Horizontal gene transfer (HGT) is the transfer of genetic material across species boundaries and has been a driving force in prokaryotic evolution. HGT involving eukaryotes appears to be much less frequent, and the functional implications of HGT in eukaryotes are poorly understood. We test the hypothesis that parasitic plants, because of their intimate feeding contacts with host plant tissues, are especially prone to horizontal gene acquisition. We sought evidence of HGTs in transcriptomes of three parasitic members of Orobanchaceae, a plant family containing species spanning the full spectrum of parasitic capabilities, plus the free-living Lindenbergia Following initial phylogenetic detection and an extensive validation procedure, 52 high-confidence horizontal transfer events were detected, often from lineages of known host plants and with an increasing number of HGT events in species with the greatest parasitic dependence. Analyses of intron sequences in putative donor and recipient lineages provide evidence for integration of genomic fragments far more often than retro-processed RNA sequences. Purifying selection predominates in functionally transferred sequences, with a small fraction of adaptively evolving sites. HGT-acquired genes are preferentially expressed in the haustorium-the organ of parasitic plants-and are strongly biased in predicted gene functions, suggesting that expression products of horizontally acquired genes are contributing to the unique adaptive feeding structure of parasitic plants.

Study Design for Sequencing Studies.

Once a biochemical method has been devised to sample RNA or DNA of interest, sequencing can be used to identify the sampled molecules with high fidelity and low bias. High-throughput sequencing has therefore become the primary data acquisition method for many genomics studies and is being used more and more to address molecular biology questions. By applying principles of statistical experimental design, sequencing experiments can be made more sensitive to the effects under study as well as more biologically sound, hence more replicable.

Selecting Superior De Novo Transcriptome Assemblies: Lessons Learned by Leveraging the Best Plant Genome.

Whereas de novo assemblies of RNA-Seq data are being published for a growing number of species across the tree of life, there are currently no broadly accepted methods for evaluating such assemblies. Here we present a detailed comparison of 99 transcriptome assemblies, generated with 6 de novo assemblers including CLC, Trinity, SOAP, Oases, ABySS and NextGENe. Controlled analyses of de novo assemblies for Arabidopsis thaliana and Oryza sativa transcriptomes provide new insights into the strengths and limitations of transcriptome assembly strategies. We find that the leading assemblers generate reassuringly accurate assemblies for the majority of transcripts. At the same time, we find a propensity for assemblers to fail to fully assemble highly expressed genes. Surprisingly, the instance of true chimeric assemblies is very low for all assemblers. Normalized libraries are reduced in highly abundant transcripts, but they also lack 1000s of low abundance transcripts. We conclude that the quality of de novo transcriptome assemblies is best assessed through consideration of a combination of metrics: 1) proportion of reads mapping to an assembly 2) recovery of conserved, widely expressed genes, 3) N50 length statistics, and 4) the total number of unigenes. We provide benchmark Illumina transcriptome data and introduce SCERNA, a broadly applicable modular protocol for de novo assembly improvement. Finally, our de novo assembly of the Arabidopsis leaf transcriptome revealed ~20 putative Arabidopsis genes lacking in the current annotation.

Detecting and Characterizing the Highly Divergent Plastid Genome of the Nonphotosynthetic Parasitic Plant Hydnora visseri (Hydnoraceae).

Plastid genomes of photosynthetic flowering plants are usually highly conserved in both structure and gene content. However, the plastomes of parasitic and mycoheterotrophic plants may be released from selective constraint due to the reduction or loss of photosynthetic ability. Here we present the greatly reduced and highly divergent, yet functional, plastome of the nonphotosynthetic holoparasite Hydnora visseri (Hydnoraceae, Piperales). The plastome is 27 kb in length, with 24 genes encoding ribosomal proteins, ribosomal RNAs, tRNAs, and a few nonbioenergetic genes, but no genes related to photosynthesis. The inverted repeat and the small single copy region are only approximately 1.5 kb, and intergenic regions have been drastically reduced. Despite extreme reduction, gene order and orientation are highly similar to the plastome of Piper cenocladum, a related photosynthetic plant in Piperales. Gene sequences in Hydnora are highly divergent and several complementary approaches using the highest possible sensitivity were required for identification and annotation of this plastome. Active transcription is detected for all of the protein-coding genes in the plastid genome, and one of two introns is appropriately spliced out of rps12 transcripts. The whole-genome shotgun read depth is 1,400× coverage for the plastome, whereas the mitochondrial genome is covered at 40× and the nuclear genome at 2×. Despite the extreme reduction of the genome and high sequence divergence, the presence of syntenic, long transcriptionally active open-reading frames with distant similarity to other plastid genomes and a high plastome stoichiometry relative to the mitochondrial and nuclear genomes suggests that the plastome remains functional in H. visseri. A four-stage model of gene reduction, including the potential for complete plastome loss, is proposed to account for the range of plastid genomes in nonphotosynthetic plants.

Comparative transcriptome analyses reveal core parasitism genes and suggest gene duplication and repurposing as sources of structural novelty.

The origin of novel traits is recognized as an important process underlying many major evolutionary radiations. We studied the genetic basis for the evolution of haustoria, the novel feeding organs of parasitic flowering plants, using comparative transcriptome sequencing in three species of Orobanchaceae. Around 180 genes are upregulated during haustorial development following host attachment in at least two species, and these are enriched in proteases, cell wall modifying enzymes, and extracellular secretion proteins. Additionally, about 100 shared genes are upregulated in response to haustorium inducing factors prior to host attachment. Collectively, we refer to these newly identified genes as putative "parasitism genes." Most of these parasitism genes are derived from gene duplications in a common ancestor of Orobanchaceae and Mimulus guttatus, a related nonparasitic plant. Additionally, the signature of relaxed purifying selection and/or adaptive evolution at specific sites was detected in many haustorial genes, and may play an important role in parasite evolution. Comparative analysis of gene expression patterns in parasitic and nonparasitic angiosperms suggests that parasitism genes are derived primarily from root and floral tissues, but with some genes co-opted from other tissues. Gene duplication, often taking place in a nonparasitic ancestor of Orobanchaceae, followed by regulatory neofunctionalization, was an important process in the origin of parasitic haustoria.